Induced decoherence and entanglement by interacting quantum spin baths 
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The reduced dynamics of a single or two qubits coupled to an interacting quantum spin bath 
modeled by a XXZ spin chain is investigated. By using the method of time-dependent density 
matrix renormalization group (t-DMRG), we go beyond the uniform coupling central spin model 
and evaluate nonperturbatively the induced decoherence and entanglement. It is shown that both 
decoherence and entanglement strongly depend on the phase of the underlying spin bath. We show 
that in general, spin baths can induce entanglement for an initially disentangled pair of qubits. 
Furthermore, when the spin bath is in the ferromagnetic phase, because qubits directly couple to the 
order parameter, the reduced dynamics shows oscillatory type behavior. On the other hand, only for 
paramagnetic and antiferromagnetic phases, initially entangled states suffer from the entanglement 
sudden death. By calculating concurrence, the finite disentanglement time is mapped out for all 
phases in the phase diagram of the spin bath. 
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I. INTRODUCTION 

Spin qubits are promising candidates for quantum in- 
formation processing due to their long decoherence and 
relaxation time. 1 ' 2 Some schemes, such as solid state spin 
qbuits, further enjoy the potential scalability via the in- 
tegration with nanotechnology. 3 However, spin qubits are 
not totally immune from the ubiquitous decoherence. To 
describe the bath that causes the decoherence of spin 
qubits, it is known that in some cases, the bath is bet- 
ter modeled by spins instead of delocalized oscillators, 
resulting in the so-called spin baths. 4 It has been argued 
that the influence of spin baths may be qualitatively dif- 
ferent from bosonic baths and non-Markovian dynam- 
ics can easily emerge.— Due to the growing interest in 
spin baths, the decoherence behavior and the entangle- 
ment dynamics of few qubits coupled to spin baths have 
been studied extensively in recent years. Early works fo- 
cus on the decoherence due to independent spins.— Here 
although the proposed model formally resembles a spin 
boson model, non-Markovian already emerges even when 
bath modes are not interacting. 7 In real baths, however, 
spins are not independent. It is therefore important to 
include effects due to interactions of spins in the bath. 
Nonetheless, the inclusion of the intra-spin interaction 
in the bath complicates the problem and only for some 
limited models with high symmetry, exact reduced dy- 
namics can be identified.* Beyond models with exact so- 
lutions, approximated dynamics was obtained by using 
mean-field^ or perturbative approaches^to handle more 
generic models. The most common model employed in 
these works is the " central spin model" , where the qubits 
are uniformly coupled to all spins of the bath. While an- 
alytical derivations are possible in these models, they are 
less realistic and are more difficult to be implemented 
experimentally. A non-perturbative approach that can 
capture the non-Markovian effects induced by interact- 
ing spin bath with generic coupling to qubits is hence 
highly desirable. 

To overcome the difficulty associated with interacting 



spins, we utilize the method of time-dependent density 
matrix renormalization group (t-DMRG)— ^ to investi- 
gate the reduced dynamics of single or two qubits coupled 
to an interacting spin chain. Recently, t-DMRG has been 
used to study the single qubit pure dephasing induced 
by a XXZ anisotropic spin chain^ The advantage of t- 
DMRG is its ability to calculate reduced dynamics even 
when the spin bath is not integrable and the coupling is 
not uniform. Due to the accumulation of errors, t-DMRG 
will eventually run away at large time 1 ^ but this does not 
impose serious limitation since for the study of quantum 
information we are mostly interested in some smaller time 
scale. In this work we apply the method of t-DMRG to 
investigate both pure dephasing and general decoherence 
model of qubits coupled to spin baths. Single qubit deco- 
herence as well as two qubits (dis)-entanglement dynam- 
ics are investigated. It is shown that both decoherence 
and entanglement strongly depend on the phase of the 
spin bath. In general, we find that spin baths can in- 
duce entanglement for an initially disentangled pair of 
qubits. However, when the spin bath is in the ferromag- 
netic phase, because qubits directly couple to the order 
parameter, the reduced dynamics shows oscillatory type 
behavior. On the other hand, only for paramagnetic and 
antiferromagnetic phases, initially entangled states suffer 
from the entanglement sudden deathj 25 i 26 

To quantify the single spin decoherence, we evaluate 
the evolution of the Loschmidt echo^ We analyze the 
relation between the short time Loschmidt echo decay 
parameter and the quantum phases of the spin bath, as 
it has been pointed out that these two are closely re- 
lated, especially when a symmetry breaking occurs in the 
bathi^ We use the temporal evolution of concurrence 1 ^ to 
study the entanglement dynamics. One important issue 
of entanglement dynamics is the possibility of creating 
entanglement through a common bath for originally dis- 
entangled qubits. It has been shown that induced en- 
tanglement via a common bath is possible for bosonic 
and fermionic baths. 17 ' 18 For spin baths such a possibil- 
ity has been explored for non-interacting spin bath 1 ^ as 
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well as interacting ones & 20 i 21 ' 22 but is restricted to uni- 
form coupling models. It will be shown later in this paper 
that induced entanglement is possible for local coupling 
model considered in this work. We note that the induced 
entanglement is also closely related to recent proposals 
of quantum communication and teleportation via spin 
chain] 23 ! 24 

Another important issue is the disentanglement dy- 
namics of an initially entangled state. It has attracted 
much attention in recent years since Yu and Eberly 25 
and Jakobczyk and Jamroz 26 predicted that two ini- 
tially entangled state without interaction can become 
completely disentangled at finite time. This feature has 
been termed entanglement sudden death (ESD). ESD has 
been studied theoretically within various models— i 28 i 29 
and has been demonstrate experimentally.— These mod- 
els, however, are restricted to Markovian bosonic bath 
or classical noise. In this work we explore for the first 
time if ESD-like phenomena can occur for spin baths in 
non-Markovian regime. In particular, we will show that 
only when the spin bath are paramagnetic or antiferro- 
magnetic, the phenomenon of ESD occurs, while when 
the spin bath is in the ferromagnetic phase, concurrence 
shows oscillatory behavior. As understanding the nature 
of the (dis)- entanglement dynamics constitutes an im- 
portant step for quantum engineering these systems, our 
results are of practical usage for future quantum infor- 
mation processing. 

The paper is organized as follows. In Sec. |TT]we present 
our model Hamiltonian and briefly discuss how to apply 
t-DMRG to analyze the model Hamiltonian. In Sec. IIIII 
we present our results of single qubit decoherence while 
in Sec. HVI the results of (dis) -entanglement dynamics are 
presented. In Sec|Vl we summarize and discuss implica- 
tion of our results. 



II. THEORETICAL FORMULATION 

We consider a system-bath model which is described 
by the total Hamiltonian H = H sys + Hbath + Hint > where 
H sys is the Hamiltonian of a single or two qubits system, 
Hbath is the Hamiltonian of a spin bath and Hi nt repre- 
sents the interaction between qubits and the bath. We 
shall set H sys = but our method can be applied to a 
generic H sys . We shall assume that the spin bath is a 
spin chain characterized by the XX Z Heisenberg model 

"Hbath = J X! ( S i S i+l + S f S i+l + &S*Si +1 ) , (1) 

where J > 0. It is known that the XXZ Heisenberg 
model has a very rich structured The system is ferro- 
magnetic for A < 1, antiferromagnetic (Ising-type) for 
A > 1, and critical (XY-type) for -1 < A < 1. It 
also encompasses the XY model where A = 0. The 
most general linear coupling between qubit A(B) and 
the bath can be expressed as H int = J2 i a efs^^Sf, 
where a = x, y, z and i = 1, . . . , N. Here characterizes 



the coupling of the spin qubit to the ith spin in the spin 
chain. For most situations, the more interesting cases are 
e < QtLlS. and hence we shall concentrate on negative e. 
Our numerical method, however, can be applied equally 
well to cases with positive e. In our work, both the Ising 
coupling and isotropic Heisenberg coupling will be con- 
sidered. The Ising coupling (ef = ef = 0) gives rise to 
a pure dephasing model while the isotropic Heisenberg 
coupling (ef = e( — ef ^ 0) induces both dephasing and 
energy relaxation. The range of the coupling is crucial 
for characterizing the interaction of qubits to the spin 
bath. For uniform coupling e; is independent of i. This is 
unrealistic but for uniform coupling the Loschmidt echo 
and the entanglement dynamics can be calculated ex- 
actly by using Jordan- Wigner transformation when the 
spin bath is the XY model (A = 0) i 20 i 32 However, the 
more realistic coupling model is the local coupling model 
in which only the coupling to the closest spin is nonvan- 
ishing. Nonetheless, there is no analytic solutions known 
for this model. The reduced dynamics is less studied 
but is more relevant to real experiments. In this case, 
if the spin bath is ferromagnetic, the qubit couples di- 
rectly to the order parameter (the magnetization); while 
if the spin bath is antiferromagnetic or paramagnetic, the 
qubit does not couple to the order parameter. Hence the 
reduced dynamics exhibits completely different behavior 
in different phases. As the local coupling model is more 
relevant to real experiments, in the following we shall 
concentrate on the local coupling model. 

We now briefly outline the procedure to evaluate the 
reduced dynamics of the qubits and other derived quanti- 
ties. For a given set of parameters, we first employ static 
DMRG 33 to find the ground state \G) of the spin chain, 
where open boundary condition (OBC) is used. We as- 
sume that at t = the initial total state is a product 
state of the form: |$(0)) = \ip sys (0))\G) , where \^ sys (0)) 
is some particular system state that we are interested in. 
Formally the evolution of the reduced density matrix can 
be obtained by evolving first the total state 

Mt))=e- tHt Hsys(0))\G), (2) 

then tracing off the spin bath 

Psys (t)=TT ba t h mt))(<P{t)\. (3) 

Loschmidt echo and concurrence then can be evaluated 
from p sys (t). In general evolving such a state is a 
formidable task. t-DMRG, however, provides a way to 
efficiently evolve such a state with high accuracy for a 
quasi-one dimensional system. We note that the degrees 
of freedom of the qubits are kept exactly during the t- 
DMRG calculation by targeting an appropriate state. 
The dimension of the truncated Hilbert space is set to 
be D — 100. For short time decay simulation we set 
JSt = 10~ 3 in the Trotter slicing while for entanglement 
dynamics we set JSt = 0.1 — 0.5 to balance the Trotter 
error and truncation error. 
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III. SINGLE QUBIT DECOHERENCE 

In this section, we present our results for the sin- 
gle qubit decoherence which is characterized by the 
Loschmidt echo. The Loschmidt echo has been used ex- 
tensively to quantify the single qubit decoherence, es- 
pecially its connection to the quantum criticality of the 
spin baths . 13 i 34 ' 35 Loschmidt echo can be understood 
intuitively as follows: Consider an initially disentan- 
gled total state (C + |+) +C_|— )) (gs |G), at some later 
time t it will evolve into an entangled state C+(t)\+) (g) 
\^f + (t))+C-(t)\-)<Si\^-{t)}. Loschmidt echo, defined as 
C{t) = \(^ + (t)\^-(t))\ 2 , clearly measure the coherence 
between |+) and |— ). When C — 1 the qubit is disentan- 
gled from the bath while when C = the qubit is totally 
entangled with the bath. 

We start by noting that for numerical calculation on 
finite length, all dynamics will show quasi-periodic be- 
havior. The quasi-period is known as the revival time. 
Since in our numerical calculation, the spin bath is a 
chain of finite length, it is essential to identify the revival 
time for each length to avoid unphysical results due to re- 
vival. As a zeroth order approximation, the revival time 
is proportional to the length and inverse proportional to 
the maximum phase velocity of the spin chain. In Figs.[T] 
and [21 we plot the Loschmidt echo as a function of time 
for the case of Ising and Heisenberg coupling using var- 
ious lengths and A, from which the revival time of the 
echo can be easily identified. We also plot the revival 
time as a function of the length. One clearly observes 
the liner dependence of revival time on the length. We 
find that for Ising coupling the minimal value Loschmidt 
echo reached is non-zero unless a very strong coupling 
strength is taken (not shown here), while for Heisenberg 
coupling, the Loschmidt echo reaches zero if the length 
is longer than some A dependent critical length. 
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FIG. 1: (a)(b)(c) C(t) as a function of time for different 
lengths and A. (d)Revival time as a function of length. The 
coupling is of Ising type with e = —0.3. 

To compute the single qubit decoherence, we couple 
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FIG. 2: (a)(b)(c) C(i) as a function of time for different 
lengths and A. (d)Revival time as a function of length. The 
coupling is of Heisenberg type with e = —0.3. 



the qubit to a single site of spin chain, which is taken 
to be the middle site of the chain to suppress bound- 
ary effects. We tune the spin bath to different quantum 
phases by changing the parameter A. In the ferromag- 
netic regime, a small uniform external field is introduced 
during the static-DMRG calculation but is turned off 
during time evolution. It is numerically checked that 
the numerical results reported below are insensitive to 
the magnitude of the applied external field. Similarly, 
when the bath is in the Ising antiferromagnetic ground 
state of the XXZ model, a small staggered external field 
is applied to lift the two-fold degeneracy in the ground 
stated 

When the spin is in the Ising antiferromagnetic regime 
or the XY critical regime, we find that in short time 
the behavior of the Loschmidt echo decay is Gaussian, 



C(t) 



where a is the decay parameter^ In Fig. [3] 



we plot the decay parameter as a function of A for the 
case of Ising coupling. In the ferromagnetic regime (A < 
— I), the qubit decay is completely suppressed (a = 0). 
This is a consequence of the Ising coupling in which both 
|+) (g) \G) and |— ) £g> \G) are eigenstates to the system and 
hence C{t) = 1. 

Clearly, the decay parameter is largest in the critical 
regime (—1 < A < I) and it decreases gradually to zero 
as one moves into the antiferromagnetic regime (A > 1). 
For single link scenario, the decay parameter is almost 
featureless within the critical regime. Our numerical re- 
sults also show that if the qubit is coupled to multiple 
sites, decay parameter acquires a weak dependence on A 
and the transition near A = 1 becomes less sharp (not 
plotted). Note that the decay parameter becomes sen- 
sitive to the magnitude of the small staggered field ap- 
plied when the spin bath is close to the phase boundary 
(A ~ I). This is due to the fact that for finite N the bar- 
rier between two degenerate ground states is finite and 
approaches zero as A approaches 1. The ground state 
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FIG. 3: Decay parameter as a function of A for the case of 
Ising coupling. Here e = —0.3, and length TV = 80. 
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FIG. 4: C{t) as a function of time when the spin bath is in the 
ferromagnetic phase for the case of Heisenberg coupling. Here 
length TV = 80 and the coupling strength e = —0.3 (upper) or 
e = 0.3 (lower) Clear oscillatory behaviors are seen for both 
positive and negative e. 



obtained by static DMRG includes a small mixture of 
the degenerate state, which is sensitive to the strength 
of the staggered field. For larger A, the barrier between 
two degenerate ground states increases as one moves deep 
into the antiferromagnetic regime. As a result, the decay 
parameter becomes less sensitive to the strength of the 
stagger field for A 3> 1. 

We now turn to the case of Heisenberg coupling. We 
first note that in the ferromagnetic regime, the qubit cou- 
ples to the order parameter. Therefore, the qubit is ef- 
fectively in an average magnetic field (Si). As a conse- 
quence of the Heisenberg coupling, the qubit will precess 
about (Si). Since magnons are generated at the same 
time when the qubit evolves, (Si) starts to deviate from 
1/2 and results in oscillations in the reduced dynamics. 
Fig. 2] shows some typical oscillating behavior of C(t) in 
this scenario. Clearly, the reduced dynamics is no longer 
Gaussian. Therefore, we shall not mark the ferromag- 
netic regime in the followings. 

In Fig. [5] we plot the decay parameter as a function 
of A for the case of Heisenberg coupling. The overall 
behavior is very similar to the case of Ising coupling ex- 
cept that the decay parameter depends weakly on A in 
the critical regime. This is different from the Ising cou- 
pling case shown above but is similar to the multiple sites 
Ising coupling case. For both Ising and Heisenberg cou- 
pling we find a discontinuity in the behavior of C(t) at 
A = — 1 and a first derivative discontinuity at A = 1. 
These discontinuities coincides with the phase boundary 
of the underlying spin chain. Different behavior at the 
A = ±1 can be traced back to the different nature of the 
ferromagnetic and antiferromagnetic transition and the 
close relation between the decoherence and the quantum 
criticality of the bath is clearly demonstrated. 



IV. ENTANGLEMENT DYNAMICS 
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FIG. 5: Decay parameter as a function of A for the case of 
Heisenberg coupling, e = —0.3, and length TV = 80 



In this section, we investigate the entanglement dy- 
namics of two qubits that couple to the spin bath. There 
are two central issues to be addressed. The first issue 
is the possibility of entanglement creation via the com- 
mon spin bath for a pair of initially disentangled qubits 
without direct interaction. The second issue is the disen- 
tanglement dynamics of an initially entangled state. In 
particular, we would like to address the issue if qubits 
influenced by spin baths also suffer from entanglement 
sudden death and if the entanglement sudden death de- 
pends on the quantum phase that qubits couple toi2ii To 
characterize the entanglement, we shall use concurrence 
as the measurement of entanglement^ For a given re- 
duced density matrix p(t), the concurrence is defined as 
C = max{Ai — A2 — A3 — A4, 0}, where Ai > A2 > A3 > A4 
are the square roots of the eigenvalues of the operator 
p(a v ® a v )p*(a v ® a v ) and p* is the complex conjugation 
of p. 
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FIG. 6: Entanglement dynamics for an initially disentangled 
pair of qubits for the case of Ising coupling Here A = and 
AT-sn 
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FIG. 7: Entanglement dynamics for an initially disentangled 
pair of qubits for the case of Heisenberg coupling. Here A = 
and N = 80. 



A. Entanglement creation 

It has been shown that entanglement can be created 
without direct interaction if two qubits interact with a 
common bosonic bathil or a fermionic bath^ Of par- 
ticular interest to us is the onset time of the entangle- 
ment, the strength of the induced entanglement, and 
the time scale where the induced decoherence eventually 
takes over. These considerations are important in deter- 
mining if such an induced entanglement is useful in real 
quantum computation. The issue is also closely related 
to the proposals of induced interaction via a common 
bath ) 38 i 39 where the effect of induced decoherence from 
the same bath is usually neglected during the derivation. 

In Fig. [6] we plot the concurrence as a function of time 
using various coupling strength and inter qubit distance. 
The coupling between qubits and the spin bath is of 
Ising type, which gives rise to a pure dephasing model. 



We assume that the coupling strength is the same for 
two qubits (ei = £2) and the initial state is taken to be 
^=(|00) + |01> + 1 10) + |11)). We shall set A = but sim- 
ilar results can ba obtained for A ^ 0. Before we discuss 
our findings in more details we would like to comment 
that if Markovian approximation or uniform coupling as- 
sumption are taken then one can no longer discuss the 
inter-qubit distance dependence. The relation between 
entanglement dynamics and inter-qubit distance, how- 
ever, is gaining interest since people begin to explore the 
non-Markovian effects of the bathi 18 i 40 i 41 We terminate 
the simulation at half of the revival time, where usu- 
ally Loschmidt echo reaches its minimum, to avoid the 
unphysical dynamics due to the revival. We also nu- 
merically check the finite size effect by comparing the 
entanglement dynamics from different chain lengths. We 
find that the results from different chain lengths agree 
with each other reasonably well. The length of the chain 
mainly set an upper bound for the simulation time. 

We find that for this configuration it is possible to 
create entanglement via the spin bath. In particular, 
for weaker coupling strength the induced entanglement 
rises more slowly but can reach a higher value; while for 
stronger coupling the induced entanglement rises more 
rapidly. The maximal concurrence reached, however, is 
lower. This is because larger coupling strength also lead 
to stronger decoherence. It is also evident from the figure 
that the entanglement creating rate decreases as the inter 
qubit distance increases, which is typical for this kind of 
induced interaction. We find that for large enough e and 
the smaller enough inter qubit distance the concurrence 
shows oscillatory behavior. In these cases, the coupling is 
strong enough to create concurrence oscillation but also 
weak enough to prevent the bath to totally disentangle 
the qubits. The delicate interplay between induced deco- 
herence and induced entanglement indicates that using 
such an induced entanglement for quantum computation 
is a tricky task. One has to tune the coupling to be within 
the right window to balance the effect from each side. 

In Fig. [7] we plot the concurrence as a function of time 
for the case of Heisenberg coupling, starting form the 
same initial condition. Qualitatively the behavior is sim- 
ilar to the case of Ising coupling. We find that the max- 
imal entanglement that can be reached is smaller. This 
is because for Heisenberg coupling the Loschmidt echo 
always decays to zero regardless the coupling strength 
while for Ising coupling the minimal Loschmidt echo 
value is a decreasing function of the coupling strength 
and is not zero. We also find that the onset time is 
roughly proportional to the inter qubit distance. This 
is expected as the excitation of spin chain, which medi- 
ate the entanglement generation, travels with finite phase 
velocity. The time the excitation reaches the other qubit 
would be proportional to the inter qubit distance. How- 
ever, the concurrence oscillation is absent, indicating that 
the induced interaction is weaker for Heisenberg coupling. 
We note that it is difficult to write down an exact form 
of the induced interaction unless Markovian approxima- 
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FIG. 8: Disentanglement dynamics for an initially disentan- 
gled pair of qubits when the bath is in the ferromagnetic 
phase. The coupling is of Heisenberg type, e = —0.3 and 
TV = 80. 

tion is taken. In general the induced interaction is time- 
dependent and is accompanied by a complicate decoher- 
cncc effect. It is, however, possible to perform quantum 
state tomography experimentally or numerically to ex- 
tract Kraus operators. The Kraus operators can then be 
used to design quantum operations without using directly 
the form of the induced interaction. 
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FIG. 9: Disentanglement dynamics for an initially disentan- 
gled pair of qubits when the bath is in the antiferromagnetic 
or the XY critical phase. The coupling is of Heisenberg type, 
e = -0.3, and TV = 80. 



B. Entanglement decay 

Here we present our results for the disentanglement 
dynamics of an initially entangled state. To investigate 
the possibility of ESD in spin bath we start from an ini- 
tial state of the form |^ sys (0)) = a|00) + /3|H), with an 



FIG. 10: Inverse finite disentanglement time as a function of 
A. Here the coupling is the Heisenberg type, e = —0.3, and 
TV = 80. 



initial concurrence C(0) = 2|a/3*|. Two qubits are set be 
20 sites apart so that the decoherence of individual qubit 
are nearly independent and the coupling is of Heisenberg 
type. We first show a typical behavior of concurrence in 
the ferromagnetic regime in Fig. [5J Clearly, as for the re- 
duced dynamics for a single qubit, the concurrence shows 
oscillatory behaviors. Hence qubits in the ferromagnetic 
regime do not suffer from ESD and the envelop of the 
entanglement decays exponentially. 

In Fig[5] we plot the disentanglement dynamics of two 
states in the XY critical and antiferromagnetic regimes. 
All of the computations start from the same initial con- 
currence with two set of coefficients, a/ ft — 1/V3 or 
ft/a = 1/ v3, corresponding to two different initial states. 
We find that in the critical regime (—1 < A < +1) 
both states suffers from ESD. Furthermore, the entangle- 
ment dynamics of two states are identical in the critical 
regime, which is due to that the rotational symmetry is 
not broken in the critical regime and 1 00} is equivalent 
to 1 11) by Zi symmetry along the quantization axis. In 
both ferromagnetic and antiferromagnetic phase, where 
the rotational symmetry is broken, we find that the en- 
tanglement dynamics for these two states starts to de- 
viate from each other. In most of the antiferromagnetic 
regime both states do not suffer from ESD. When A is 
close to the phase boundary, however, both of two states 
that correspond to a/ ft — l/y/3 and ft/a = l/\/3 suf- 
fer from ESD and have slightly different disentanglement 
time. In Fig. [TU] we plot the inverse of disentanglement 
time, which is defined as the time when concurrence be- 
comes zero, as a function of A. Starting from A = — 1, it 
shows monotonic decrease. Across the phase boundary, 
A = 1, the inverse of disentanglement time develops a 
small bump and persists into the antiferromagentic phase 
and finally decreases to zero around A = 1.2. The exis- 
tence of finite region with finite disentanglement time in 
the antiferromagnetic regime is due to that when A ap- 
proaches 1, the barrier between two degenerate ground 
states approaches zero. For finite TV and finite value of 
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A — 1, the ground state obtained by static DMRG in- 
cludes a small mixture of the degenerate state so that it 
resembles the XY critical state and results in finite disen- 
tanglement time. In the thermodynamic limit (N — > oo), 
the region with finite disentanglement time in the antifcr- 
romagnetic regime shrinks down, resulting discontinuity 
at the phase boundary A = 1. The overall behavior 
found in the above is different from those reported in 
the early work22 in which two states investigated in the 
above posses different entanglement dynamics and only 
one of them suffer from ESD. The difference is due to 
that the model adopted in Rcf. 27 includes the effect of 
spontaneous decay, breaking the symmetry between |0) 
and 1 1); while for spin bath in XY critical regime such 
symmetry breaking is absent. It is important to note that 
there exits a subregime, roughly around —0.3 < A < 0.6, 
in which entanglement shows revival after some dark pe- 
riod. Note that the entanglement revival after some dark 
period was also reported in Ref42 where a photonic mul- 
timode vacuum bath is assumed and the revival is at- 
tributed to the two photon decay. In our work, the origin 
of the revival is less clear. We believe that the existence 
of such sub-regime is due to the competition between en- 
tanglement decay and Loschmidt echo decay. Within this 
subregime the Loschmidt echo decay more slowly, giving 
the system a chance to revive after the first ESD. 

V. CONCLUSION 

In summary, the decoherence and (dis)-entanglement 
dynamics induced by spin baths are investigated non- 
perturbatively by using t-DMRG. For both pure de- 



phasing model (Ising coupling) and general decoherence 
model (Heisenberg coupling) we calculate the short time 
decay parameter of the Loschmidt echo. We find that in 
both cases the decay parameter is closely related to the 
phase of the underlying spin chain. In the ferromagnetic 
regime, the reduced dynamics shows oscillatory behavior 
while in the XY critical and antiferromagnetic regimes, 
the decay parameter shows a first derivative discontinu- 
ity at A = +1. We evaluate the entanglement dynamics 
of a pair of initially disentangled qubits which are close 
to each other. We demonstrate that it is possible to in- 
duced entanglement via their common interaction with 
the spin bath. The competition between induced deco- 
herence and entanglement can be easily seen in the cou- 
pling strength dependent behavior of the entanglement 
onset time, growth rate, and the maximal entanglement 
reached. Finally we investigate the disentanglement dy- 
namics of a pair initially entangled qubits which are far 
from each other. For the two initial states we studied, 
we find that their disentanglement dynamics are identi- 
cal and suffer from ESD in the critical regime. Their 
disentanglement dynamics begin to deviate from each 
other in both the ferromagnetic and the antiferromag- 
netic regime. They no longer suffer from ESD in the fer- 
romagnetic regime but still suffer from ESD if the chain 
is near the antiferromagnetic transition. It is shown that 
the inverse of finite disentanglement time has a close rela- 
tion to the phase of the spin bath and shows monotonic 
decrease behavior as one moves into the antiferromag- 
netic regime. 

We thank Prof. Hsiu-Hau Lin for inspiring discussions 
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ence Council of Taiwan. 
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